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SUMMARY 

,/ 

This  is  the  final  technical  report  on  a study  of  finite  amplitude  attenua- 
tion ^hat  has  been  performed  under  Contract  N00039-75-C-0263,  Program  Code 
5G10,  ARPA  Order  2910.  The  program  objective  has  been  to  determine  by 
experimental  and  theoretical  means  the  limitations  imposed  by  finite  amplitude 
attenuation  on  the  performance  of  parametric  sonars.  The  results  are  intended 
for  use  in  assessing  the  potential  of  parametric  receivers  and  transmitters  in 
military  applications. 

The  axial  response  of  the  parametric  receiving  array  has  been  calculated 
in  two  ways.  The  first  is  by  use  of  a one  dimensional  model  in  which  the 
secondary  signal  levels  are  predicted  from  weak  shock  theory.  The  second 
is  Dy  numerical  evaluation  of  the  scattering  integral  solution  in  which  finite 
amplitude  attenuation  is  introduced  by  means  of  a taper  function  derived  from 
the  one  dimensional  solution. 

An  experiment  has  been  performed  in  which  the  axial  secondary  signal 
level  as  a function  of  pump  source  level  was  measured.  These  data  agree  with 
the  theoretical  results,  but  the  agreement  is  not  conclusive  because  the  maxi- 
mum pump  source  level  was  not  sufficiently  high.  However,  a comparison 
with  data  from  another  investigation  confirms  the  theory  for  cases  when  the 
receiving  hydrophone  is  well  in  the  farfield  of  the  pump. 

Comparison  with  a set  of  experimental  results  for  which  the  interaction 
was  almost  entirely  in  the  pump  nearfield  indicates  that,  unlike  the  farfield 
case,  the  nearfield  behavior  may  be  quite  different  from  that  predicted  by  the 
present  theory.  Caution  is  therefore  advised  in  applying  the  present  results 
to  nearfield  interactions. 
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In  addition  to  the  parametric  receiving  array  work,  several  parametric 
transmitting  array  taper  functions  have  been  compared  in  the  context  of 
numerical  evaluations  of  the  scattering  integral.  A mean  taper  function  is 
derived  from  weak  shock  theory  and  is  compared  with  the  others. 

In  all  cases  the  results  are  surprisingly  similar  to  those  obtained 
using  the  model  g'ven  by  Mellen  and  Moffett. 
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1.  o INTRODUCTION 

This  is  the  Final  Technical  Report  under  Contract  N00039-75-C- 
02  63,  Program  Code  5G10,  ARPA  Order  2910.  The  program  objective  has 
been  to  determine,  by  experimental  and  theoretical  means,  the  limitations 
imposed  by  finite  amplitude  attenuation  on  the  performance  of  parametric 
sonar.  The  main  emphasis  in  this  work  has  been  on  the  parametric  receiving 
array  which,  in  contrast  to  the  parametric  transmitting  array,  has  received 
only  a small  amount  of  previous  attention  so  far  as  finite  amplitude  attenua- 
tion is  concerned. 

The  axial  response  of  the  parametric  receiving  array  has  been 
modeled  in  two  ways.  The  first  model  consists  of  a one  dimensional  wave 
undergoing  time  domain  distortion  to  produce  sum  and  difference  frequencies 
as  it  propagates  from  the  pump  projector  to  the  receiving  hydrophone.  In 
the  second  model,  the  parametric  receiver  is  represented  as  a volume  source 
of  sum  and  difference  frequencies  that  scatter  from  each  point  in  the  inter- 
action volume  and  are  summed  at  the  receiving  hydrophone.  In  the  latter 
model  the  finite  amplitude  attenuation  is  introduced  by  means  of  a taper 
function  derived  from  the  one  dimensional  model. 

With  the  introduction  of  a simple  correction  to  the  one  dimen- 
sional model  to  account  for  nearfield  diffraction  effects,  the  two  models 
agree  quite  well  with  each  other.  The  experimental  results  for  cases  when 
the  receiving  hydrophone  is  well  in  the  farfield  of  the  pump  are  also  in  close 
agreement  with  the  theory.  There  is,  however,  evidence  that  the  taper  model 
used  may  not  be  adequate  for  describing  the  effects  of  finite  amplitude  attenua- 
tion in  the  nearfield  of  the  pump.  There,  the  loss  due  to  saturation  (based 
on  the  only  available  set  of  experimental  data)  appears  to  be  considerably 
less  than  expected.  This  agrees  with  the  results  of  an  alternative  taper 
function  that  is  simply  equal  to  the  fundamental  component  amplitude  in 
the  saturating  pump  signal. 

The  taper  function  derived  herein  is  based  on  farfield  considera- 
tions and  is  believed  to  be  realistic.  While  it  is  tempting,  based  on  the 
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apparently  different  behavior  in  the  nearfield  to  combine  the  two  types  of 
taper,  this  has  not  been  done  owing  to  the  lack  of  theoretical  justification 
and  the  availability  of  only  a single  source  of  nearfield  data.  Thus,  appli- 
cation of  the  theory  presented  here  for  the  axial  response  of  the  parametric 
receiving  array  should  be  restricted  to  cases  in  which  the  receiving  hydro- 
phone is  well  in  the  farfield  of  the  pump. 

The  experimental  portion  of  this  investigation  produced  data  for 
the  axial  response  of  the  parametric  receiving  array  that  agrees  well  with 
the  theory  as  calculated  by  the  methods  described  herein.  However,  efforts 
to  produce  a pump  signal  of  sufficient  intensity  to  make  the  results  regard- 
ing the  effects  of  saturation  conclusive  were  unsuccessful.  It  has,  therefore, 
been  necessary  to  rely  on  results  obtained  by  other  investigators  for  valida- 
tion of  the  theory. 

As  a secondary  task  of  the  present  investigation  the  parametric 
transmitting  array  was  modeled  by  use  of  numerical  evaluation  of  the  volume 
scattering  integral  with  inclusion  of  a taper  function  capable  of  accounting 
for  finite  amplitude  attenuation  as  a function  of  source  amplitude,  range  and 
angle  relative  to  the  acoustic  axis  of  the  projector.  The  model  was  so  con- 
structed that  a variety  of  taper  functions  could  be  used  for  comparison 
purposes  in  addition  to  the  one  derived  in  this  report.  Several  such  taper 
functions  are  compared  both  as  to  their  mathematical  form  and  as  to  their 
performance  in  the  scattering  integral  evaluation.  The  performance  of  the 
model  with  the  taper  function  derived  in  this  report  is  tested  by  comparison 
of  computed  results  with  experimental  data  obtained  by  Mellenbruch  and 
Muir1  for  the  case  of  a parametric  array  reflected  from  a pressure  release 
surface  a short  distance  from  the  projector. 

The  technical  body  of  this  report  is  organized  into  two  major 
sections  following  the  introductory  material  in  Section  1.0.  Section  2.  0 
contains  the  treatment  of  the  parametric  receiving  array  including  theory, 
experiment,  comparison  with  available  data  and  conclusions.  Section  3.  0 
deals  with  finite  amplitude  effects  in  parametric  transmitting  arrays. 
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The  hnical  body  is  now  preceded  by  a brief  review  of  the 
background  of  th  ''ametric  array. 

1.  1 Background 

The  parametric  array  had  its  beginnings  about  1960  with  the 
theoretical  work  of  Westervelt.  ’ Westervelt's  analysis  showed  that  if  a 
sound  wave  consisting  of  two  discrete  high  frequencies  could  be  confined 
to  a narrow,  collimated  beam,  the  two  components  would  interact  to  produce 
a sound  wave  at  the  difference  frequency.  This  is  the  basis  of  the  parametric 
transmitting  array.  He  showed  that  the  directivity  of  the  difference  frequency 
sound  wave  far  from  the  interaction  region,  would  depend  only  on  the  attenua- 
tion at  the  primary  frequencies  and  on  the  difference  frequency  wave  number. 
The  interaction  region  is  limited  by  the  attenuation  at  the  primary  frequencies. 

In  addition  to  the  assumption  that  the  primary  beam  acted  as  a 
line  source,  Westervelt  assumed  that  the  primary  wave  attenuation  was 
adequately  described  by  linear  theory  (i.  e. , that  the  amplitude  was  not  so 
high  that  finite  amplitude  attenuation  became  a factor).  Early  experimental^ 
results  confirming  Westervelt's  theory  were  published  by  Beilin  and  Beyer. 
Westervelt  also  observed  that  a higii  frequency  pump  could  be  made  to  function 
as  a receiver  by  nonlinear  interaction  with  a low  frequency  signal.  This  is 
the  basis  of  the  parametric  receiving  array. 

Numerous  authors  have  extended  the  basic  parametric  trans- 
mitting and  receiving  array  theories  to  describe  more  realistic  geometries 
and  to  account  for  high  amplitudes.  There  have  been  two  basic  approaches 
to  handling  geometry.  The  first  is  to  approximate  the  interacting  signals 
as  one  dimensional  propagating  waves  (see,  for  example.  Refs.  5 and  6). 

The  second  is  to  perform  a three  dimensional  integration  as  was  done  by 
Muir  and  Willette.  7 The  latter  approach  is  the  more  general  but  is  less 
amenable  to  extension  to  account  for  finite  amplitude  effects.  It  has,  in 
fact,  been  the  former  approach  that  has  been  used  almost  exclusively  in 
dealing  with  finite  amplitude  effects. 
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Very  little  work  has  been  done  in  analyzing  finite  amplitude 
etiects  in  parametric  receiving  arrays.  Experimental  data  showing  satura- 
tion in  parametric  receivers  have  been  published  by  Berktay  and  Al-Temimi 
and  by  Konrad,  Mellen  and  Moffett.  A theoretical  model  was  presented 
by  Bartram  ^ for  cases  when  both  pump  and  signal  are  plane  waves.  Other 
geometries  have  not  been  dealt  with. 

Other  work  related  to  the  parametric  receiving  array  has  been 
done  by  Fenlon^in  addressing  a multiple  frequency  parametric  source. 
Schaffer  and  Blackstock  also  performed  related  work  in  their  study  of 
modulation  of  an  intense  high  frequency  signal  by  a low  frequency  signal. 


Finite  amplitude  effects  in  parametric  transmitting  arrays  have 

13  14 

been  dealt  with  by  Mellen  and  Moffett,  Merklinger,  Bartram  and 
Westervelt, 15  and  Fenlon.16’  17  In  all  cases,  quasi  one  dimensional  models 
were  utilized  although  Merklinger  did  formulate  the  problem  in  more  generality. 

The  treatment  of  finite  amplitude  attenuation  in  parametric  arrays 

is  generally  based  on  one  dimensional  models  and  often  on  theories  originaJ.lv 

formulated  for  single  frequency  sources.  Useful  theoretical  treatments  of 

18  19 

single  frequency  sources  have  been  given  by  Blackstock.  ’ A numerical 

20 

solution  for  spherical  waves  of  arbitrary  spectrum  has  been  given  by  Fenlon. 

Single  frequency  waves  from  piston  sources  have  been  studied  by  Merklinger, 

21  22  2324 

Berktay  and  Safar,  Shooter,  Muir  and  Blackstock,  Lockwood,  ’ Lockwood, 

2 r 26 

Muir  and  Blackstock,  ^ and  Browning  and  Mellen. 

Experimental  data  that  demonstrate  the  performance  of  various 

7 

experimental  parametric  transmitters  have  been  reported  by  Muir  and  Willette, 
Muir  and  Blue,^7  Merklinger,  ^ and  numerous  others.  One  specific  experi- 
ment that  is  quite  valuable  for  verifying  the  detailed  behavior  of  a parametric 
transmitting  array  theory  is  the  phase  reversed  parametric  array  experiment 
reported  by  Mellenbruck  and  Muir.1  There  a parametric  array  was  reilected 
from  a pressure  release  surface  at  a range  of  four  yards.  The  reflection 
caused  a return  of  energy  from  the  harmonics  to  the  fundamental  component 
of  the  carrier  signal,  thereby  supposedly  prolonging  the  interaction  and  in- 
creasing the  efficiency.  The  expected  efficiency  increase  was  never  measured. 
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A theoretical  attempt  by  Lockwood  at  describing  the  phase  reversed  array 
disclosed  a weakness  in  the  concept  on  which  the  model  used  was  based 
and  led  to  the  development  of  the  taper  function  model  described  in  Section 
3.  1 of  the  present  report. 


2.0 


PARAMETRIC  RECEIVING  ARRAY  AXIAL  RESPONSE 


The  prediction  of  fini  e amplitude  behavioi  in  the  parametric 
receiving  array  is,  in  the  present  work,  based  on  analysis  of  the  time 
domain  distortion  of  a composite  wave  that  initially  consists  of  an  intense, 
monochromatic  pump  wave  and  a low  amplitude  signal  wave.  As  the  com- 
posite wave  travels,  distortion  components  are  produced  by  the  modulation 
of  the  pump  wave  by  the  signal  wave  and  of  the  signal  wave  by  the  pump 
wave,  the  latter  being  generally  a small  effect. 

The  distortion  components  of  interest,  i.  e. , sum  and  difference 
frequencies,  are  generally  close  in  frequency  to  the  pump  signal,  and  in 
the  pump  farfield  they  travel  with  the  pump  wave.  On  the  other  hand,  dis- 
tortion signals  produced  in  the  pump  nearfield  are  subject  to  diffraction. 

It  is  because  of  these  nearfield  generated  signals  that  a purely  one  dimen- 
sional model  does  not,  in  general,  adequately  describe  the  axial  response 
of  the  parametric  receiver. 


The  geometric  spreading  of  the  secondary  signals  is  handled 
by  two  different  methods.  First,  the  one  dimensional  model  derived  in 
Section  2.  1 is  corrected,  as  described  in  Section  2.1.1,  by  a simple  geo- 
metrical factor  that  approximately  compensates  for  diffraction  in  the  near- 
field. Then  in  Section  2.2  a second  model  is  described  wherein  the  geometry 
is  dealt  with  rigorously  using  the  scattering  integral  approach.  Here  the 
finite  amplitude  effects  are  handled  consistent  with  the  one  dimensional 
model  by  use  of  a taper  function  derived  in  Section  2.2.  1 by  differentiating 

the  one  dimensional  solution. 

The  results  are  compared  with  experimental  data  of  other  investi- 
gators in  Section  2.  3.  Then  in  Section  2.4,  the  experimental  work  performed 
in  the  present  investigation  is  reported.  Finally,  the  conclusions  presented 
in  Section  2.  5 bring  the  discussion  of  the  parametric  receiving  array  to  a 

close. 
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2.  l Weak  Shock  Theory  Formulation  of  Parametric 

Receiving  Array  Axial  Response 

29 

In  the  first  quarterly  progress  report  of  this  investigation  a den 
vation  was  given  of  a math  model  f parametric  receiving  array  axial  response. 
The  signal  was  assumed  to  be  a plane  wave  and  the  pump  was  assumed  to  be 
a one  dimensional  pxcve,  cylindrical  or  spherical  wave.  Of  primary  interest 
are  the  plane  and  spherical  wave  formulations. 


The  derivation,  as  given,  suffers  from  two  problems.  First,  the 
secondary  frequencies  are  shown  to  be  produced  by  modulation  of  the  pump 
by  the  signal  and  also  by  modulation  of  the  signal  by  the  pump  but  in  the 
latter  case  the  secondary  frequencies  produced  appear  to  propagate  with  the 
signal.  Actually,  when  the  pump  is  spherical  iT  modulates  the  signal  along 
a spherical  wave  front  and  the  interaction  products  should  spread  spherically, 
not  as  the  plane  wave  signal.  This  fact,  which  is  apparent  from  the  scatter- 
ing integral  formulation,  means  that,  even  in  this  simple  geometry,  deter 
mination  of  the  parametric  receiving  array  axial  response  is  not  a one  dimen- 
sional problem.  The  second  problem  is  that  in  attempting  to  account  for 
absorption  by  using  modified  weak  shock  theory,  the  absorption  at  the 
secondary  frequencies  cannot  be  correctly  represented  and  must  be  introduced 
at  the  end  of  the  derivation  in  ad  hoc  fashion. 

The  derivation  as  now  presented  is  made  rigorous  by  stipulating 
that  the  signal  frequency  be  much  smaller  than  the  pump  frequency,  in  which 
case  the  sum  and  difference  frequencies  are  approximately  equal  to  the  pump 
frequency.  The  result  of  modulation  of  the  signal  by  the  pump  may  now  be 
ignored.  Two  additional  changes  are  introduced  in  the  present  derivation. 

The  signal  is  assumed  to  spread  with  a large  but  finite  radius  of  curvature 
and  the  signal  is  assumed  to  suffer  no  significant  loss  due  to  absorption 
within  the  region  of  interaction. 

With  the  stated  assumptions,  the  pump  and  signal  wave  field 

is  given  by 
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represents  the  total  retardation  of  the  phase  of  a given  wavelet  resulting  from 
the  integrated  effect  of  the  disparity  between  the  phase  velocity  and  the  small 
signal  sound  speed  (group  velocity). 

Equation  (1)  may  be  rewritten  with  the  substitutions 
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as  V = sin  0 + 6 (r)  sin  fi  0 . 

The  formulation  given  above  departs  from  the  Earns  haw  solution 
for  the  phase  distortion  of  finite  amplitude  waves  by  the  inclusion  of  the 
absorption  term.  However,  when  absorption  is  included  in  this  manner  the 
continuous  portions  of  a wave  containing  shocks  may  be  described  as  in 
weak  shock  theory  (Biackstock,  Ref.  18)  with  the  result  that  the  high  and 
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low  intensity  results  are  matched  in  a single  solution.  Now, 
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one  obtains  the  expression  for  0 as 


0 = y + a sin  0 + a sin  ft  0 

s 

(Based  on  reasoning  given  in  Ref.  2 4,  it  may  be  that  for  long  receiver  ranges 
chq  in  Eq.  9 should  be  replaced  by  2 . Normally,  however,  the  difference 
is  very  small  and  for  the  present  the  form  given  is  satisfactory.) 

It  is  now  observed  from  Eq.  (11)  that 

j 0 - y < a + a 
I 1 l - s 

Now,  if  ft  (a  + a ) < 1,  the  maximum  phase  error  incurred  by  replacing  0 by 
s 


y in  the  last  term  of  Eq.  (11)  is  1 radian.  Thus,  the  substitution  can  be 
made  without  appreciable  effect  on  the  solution  and  the  interaction  com- 
ponents may  be  obtained  from  Eq.  (7)  with  0 given  by 


0 = y + a sin  0 + a sin  ft  y 

b 


(13) 


With  the  substitution 

y*  = y + ag  sin  ft  y 


Eq.  (13)  becomes 


0 = y*  + a sin  0 


(14) 

(15) 


Then  the  first  term  of  Eq.  (7)  assumes  the  form  of  the  Eamshaw  solution. 

1 8 

Weak  shock  theory  as  formulated  by  Blackstock  may  now  be  applied.  The 
solution  expressed  as  a Fourier  series  is 


00 


sin  0 


(a)  sin  n y' 


n = 1 


(16) 
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where  the  values  for  B (o)  are  those  derived  by  Blackstock.  Eqs.  (15) 

n 

and  (16)  may  be  substituted  into  Eq.  (7)  to  give  an  explicit  solution  in  terms 
of  y*  as 


00 


V = 


■ L 

n = 1 


Bn  (a) 


sin  n y*  + 6 (r)  sin  ft 


y*  + a 


oo 

V 

n = 1 


Bn(a)  sin  n y'| 


(17) 


It  is  now  necessary  only  to  extract  the  sum  and  difference  fre- 
quency components  which  will  be  the  coefficients  of  sin  (1  ± ft)  y.  With 
the  substitution  of  Eq.  (14)  into  Eq.  (17)  and  with  the  assumption  that  ag 
is  small  and  further,  that  5 (r)  ft  a is  negligible,  it  is  possible  to  write 

b 


00 


V = ) B (a)  sin  n (y+a  sin  fty)  + 6 (r)  sin 

n s 


ft 


n = 1 


00 


y+a  I Bn(o)  sin  n y 
n = 1 


(18) 
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The  only  terms  of  the  sums  in  Eq.  (18)  that  significantly  affect 
the  sum  and  difference  frequencies  are  the  n = 1 terms.  Thus,  discarding 
the  rest  we  have 


V = B ,(cr)  sin  (y  + a sin  ft  y)  + 6(r)  sin  ft  (y+a  B^a)  sin  y)  . 
1 s 


(19) 


Equation  (19)  may  be  thought  of  as  the  sum  of  two  broadband  FM  signals. 
The  first  term  represents  the  modulation  of  the  pump  by  the  signal  and  the 
second  represents  the  modulation  of  the  signal  by  the  pump. 

Equation  (19)  may  be  written  in  exponential  notation  as 


V = Re  < 


_iy  -ia  sin  ft  y -i^y  -ifi  aB^(a)siny 

iBjWe  e S + 6(r)  i e e 


l . (20) 


Then,  using  the  Fourier  Bessel  expansion  we  obtain 


00 


V = Re  < 


. v-1 

iBAo)e~1V  > 


-imfiy 


m s 


m=  -oo 


oo 


+ 6(r)ie'iny  ) I, 
k=  -oo 


£2  a B^(a) 


-iky 


Selecting  only  the  m - ± 1 and  k = ± 1 terms  that  give  sum  and 
difference  frequencies  and  returning  to  real  notation,  we  obtain 


V = Bx(a)  J1  (ag) 


sin  ( 1 + f2)  y — sin  (1  - ^)  T! 


,1 


+ 6(r)  J,  n a B.(ct)  ' , sin  (1  + ^)  y + sin  ( 1 - «)  y 

* L 1 j i. 


(2  1) 


(22) 
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It  is  now  apparent  that  if  ft  « 1,  the  second  term  above  may  be  neglected 
in  comparison  to  the  first.  Thus,  the  result  of  modulation  of  the  pump  by 
the  signal  is  retained  and  the  result  of  the  modulation  of  the  signal  by  the 
pump,  with  its  geometrical  inconsistency's  discarded.  The  result  for  the 
sum  and  difference  frequencies  becomes 


V±  = B^o)  J1  (ag)  sin  ( 1 ± SI)  y . (23) 

12 

A similar  result  was  obtained  by  Schaffer  and  Blackstock  for  the  case  of  an 
intense  low  frequency  modulating  a lc  > amplitude  high  frequency.  However, 
in  the  present  case  it  has  been  possible  to  include  the  effects  of  shocks  in 
the  pump  wave  through  the  Bj  term  derived  from  weak  shock  theory. 

In  dimensional  notation  the  magnitudes  of  the  sum  and  difference 
pressure  signals  become 


1 _ 1 

2 

r 

o 

p, 

= P c e 

1 ±1 

o o 

l r 

n 


Bj(o)  Jj(os) 


(24) 


Normally,  a is  small  so  that  the  sma’l  argument  approximation  of  the  Bessel 
s 

function  is  valid  and  one  may  substitute 


fie  k 5 

(l  + ri> 


L + r 
IT7o 


(25) 


If  the  reference  distance  r^  is  taken  to  be  the  receiver  location,  then  r - r^  . 

Furthermore,  if  L is  large  the  value  of  o reduces  to  the  plane  wave  signal 

s 

result 


*3  ~ 60(r-r0) 

fie  k 5 r 
o o 


The  low  amplitude  limit  of  Eq.  (24)  is  then 


(26) 
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wh  ch  agrees  with  the  results  of  Berktay  and  Shooter  and  Berktay  and 
Al-Temimi3  for  the  cases  n = 1 and  n = 0,  respectively,  except  that  here 
u and  a ± have  been  approximated  by  and  <*Q  . 

2.  l.  l Use  of  the  Weak  Shock  Theory  Results  to  Calculate  Parametric 

Receiver  Amplitude  Response  - The  expressions  derived  in  the  previous 
section  may  be  used  to  calculate  the  amplitude  response  of  parametric  receiv- 
ing arrays  that  have  interaction  volumes  either  entirely  in  the  pump  nearfield 
or  predominantly  in  the  pump  farfield.  For  intermediate  cases  where  the 
receiver  is  in  the  farfield  of  the  pump  but  there  is  a significant  contribution 
from  the  nearfield,  the  plane  wave  and  spherical  wave  solutions  may  be 
combined.  However,  it  is  necessary  to  correct  the  nearfield  portion  to 
account  for  spreading  of  the  secondary  signals.  This  is  done  by  adding  to 
the  spherical  wave  solution  the  plane  wave  solution  multiplied  by  the  factor 


F = 1 \ f 

R 0 


dx  + 


where  the  quantities  (R-  RQ)  are  replaced  by  RQ  in  the  event  that  R RQ  R0  • 
Similar  factors  were  derived  in  Ref.  29  for  the  nearfield  of  the  parametric 
transmitting  array. 

The  two  terms  of  the  factor  represent  the  portion  of  the  nearfield 
beyond  one  Rayleigh  distance  from  the  receiver,  the  contributions  from  which 
undergo  spherical  spreading,  and  the  portions  within  one  Rayleigh  distance 
which  do  not.  The  initial  factor  l/R  cancels  the  range  dependence  of  the 
plane  wave  solution  which  is,  in  effect,  replaced  by  the  quantity  in  brackets. 


> • ••• 
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The  farfield  contribution  is  calculated  by  using  the  spherical 
wave  solution  with  a source  radius  of  RQ  and  an  initial  distortion  based  on 
the  value  of  a calculated  at  RQ  for  the  nearfield  distortion. 

2.  2 Scattering  Integral  Evaluation  of  Parametric 

Receiving  Array  Axial  Response 

Because  the  determination  of  the  secondary  signal  amplitude  in  a 
parametric  receiving  array  with  a piston  pump  transducer  is  not  generally  a 
problem  involving  only  one  spatial  dimension,  it  is  often  desirable  to  use 
the  scattering  integral  approach.  In  order  to  account  for  saturation  effects 
it  is  then  necessary  to  insert  a taper  function  into  the  integral.  Such  a 
taper  function  is  now  derived. 

2.2.1  Taper  Function  Derivation  - It  is  tempting  to  use  the  pump  funda- 

mental harmonic  coefficient  B^(ct)  as  the  taper  function  because  this  would 
make  the  source  strength  density  proportional  to  the  local  pump  amplitude. 
However,  that  approach  would  not  include  the  diminis hnent  of  the  secondary 
signal  due  to  pump  saturation  as  seen  in  the  v/eak  shock  solution.  There- 
fore, the  alternative  approach  of  deriving  the  taper  function  from  the  weak 
shock  result  is  adopted. 

The  finite  amplitude  taper  is  calculated  from  the  weak  shock 
solution  by  taking  the  range  derivative  of  the  high  intensity  solution  to  obtain 
the  one  dimensional  source  strength  density  and  then  dividing  by  the  low 
amplitude  source  strength  density  to  normalize  the  result.  Since  only  the 
high  intensity  taper  is  desired  the  spherical  wave  solution  is  range  normalized 
before  differentiating. 

The  two  expressions  (neglecting  absorption)  are,  from  the  pre- 
vious section  with  ot  = a - 0,  for  the  plane  wave  region  of  the  pump, 

o s 

p±  = poco2  £ Bl(ff)  tlK1  • 


and  for  the  spherical  wave  region  of  the  pump 
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t~  p±  = p0  co2  e 3l(rf  hi°J 


where 


and  a is  given  by 


a ~ /3  e k 6 
s o o 


a =/3  e k r 
o 


in  the  plane  wave  region,  and  by 


(30) 


(31) 


(32) 


a = Pe  kQ  Rq  fl  + in  (r/RQ)  ] 


(33) 


I] 


in  the  spherical  wave  region. 

By  differentiating  the  above  expressions  for  P±  and  (r/r()Pi) 
with  respect  to  r,  and  normalizing  the  plane  and  spherical  wave  results 
by  the  respective  low  intensity  limits,  the  two  taper  functions  are  obtained 
and  are  written, 

d (a  Bx) 

T = j 

p da 

for  a plane  wave,  and 

TS  = Bi  + Bi (<7>  °o 

for  a spherical  wave. 

The  above  taper  functions  are  used  in  the  numerical  evaluation 
of  the  scattering  integral  as  discussed  in  the  next  section.  They  differ 
from  the  simpler  approach  of  using  B^(o)  as  the  taper  function  in  that  they 
include  the  "absorption"  of  the  secondary  signal  by  the  saturating  pump. 

It  is  noted  however  that  the  tapers  were  derived  under  the 
assumption  that  the  pump  and  the  secondary  signals  were  propagating  in  a 
one  dimensional  fashion.  In  the  case  of  a nearfield  generated  secondary 


(34) 


(35) 
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signal,  which  undergoes  diffraction  in  passing  to  the  farfield,  it  is  not  clear 
that  the  above  tapers  are  valid.  This  question  will  be  discussed  further  in 
the  light  of  published  experimental  results  in  Section  2.3. 

2.2.2  Numerical  Evaluations  of  the  Scattering  Integral  - The  use  of 

numerical  volume  integration  in  combination  with  the  finite  amplitude  taper 
function  is  now  described.  The  programming  stems  from  a modification  of  the 
parametric  transmitting  array  program  utilized  in  this  study  which  is,  in  turn, 
a modification  of  the  Muir-Willette 7 program.  Because  this  investigation  deals 
with  describing  axial  response  no  attempt  has  been  made  to  include  oif-ax_s 
response  in  the  programs.  However,  the  inclusion  of  the  taper  function  in 
a program  that  handles  off-axis  geometries  would  be  quite  straightforward. 

Actually,  two  different  programs  have  been  developed  to  evaluate 
the  scattering  integral.  In  the  first,  a volume  integration  over  two  coordinates 
is  performed.  In  the  second,  closed  form,  high  wave  number  approximations 
are  used  to  reduce  the  solution  to  a single  integration  along  the  axis  of  the 
pump. 


The  geometry  on  which  the  programs  are  based  is  shown  in  Figure 
1.  It  consists  of  a spherical  wave  signal  source  at  an  arbitrary  distance 
(usually  long)  from  the  receiving  hydrophone  and  a pump,  assumed  for  con- 
venience in  programming  to  be  a circular  array  with  a plane  wave  nearfield 
and  a spherical  wave  farfield.  The  nearfield  of  the  pump  projector  is  taken 
to  be  a collimated  plane  wave  with  an  abrupt  transition  to  spherical  wave 
at  the  Rayleigh  distance  (area/X.)  and  an  accompanying  90  phase  shift. 

The  volume  is  formulated  in  cylindrical  coordinates  in  the  near- 
field and  in  spherical  coordinates  in  the  farfield.  The  integral  expression 
for  the  axial  secondary  pressure  is 
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(36) 


d0  dr' 


where  T is  the  nearfield  taper  and  T is  the  farfield  taper  given  by  Eqs.  (34) 
P s 

and  (35)  respectively.  The  signal  projector  to  field  point  ranges  are 


R . = V(  L + r')2  + x2  , 
si 

and 

R = /l2  + r'2  + 2L  r'  cos  0 , 

s 

while 

Rx  = V x2  + (r  - r')2 

is  the  source  to  field  point  range  in  the  nearfield,  and 


R = V r2  + r'2  - 2r  r*  cos  0 


(37) 


(38) 


(39) 


(40) 


is  the  source  to  field  point  range  in  the  farfield.  The  lower  limit  in  the  range 

integration  is  r . , the  shadow  length  of  the  pump  transducer  given  by 
3 min 

r . = area/\  . 

min  s 

The  basic  parametric  receiver  program  utilizes  a two  dimensional 
integration  of  the  above  equation.  Unfortunately,  the  phase  variations  of  the 
integrand  are  so  rapid  at  high  pump  frequencies  that  it  is  extremely  difficult 
and  expensive  to  obtain  convergence  of  the  integral. 
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To  overcome  this  difficulty  a high  wave  number  approximation  of 
the  0 integration  and  an  approximate  closed  form  x integration  have  been 
incorporated  into  a second  version  of  the  program.  This  single  integral  program 
has  been  found  to  give  equivalent  results  to  those  of  the  double  integration 
but  with  considerably  improved  convergence. 

The  angle  integration  is  approximated  as  follows:  let 


Koo 


2 (kQ  a sin  0) 

k a sin  0 
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-a  R ± ik  R 
± s s 


R 


sin  0 d0 


(41) 


Provided  that 


2 , ' ,1 

+ 2Lr  cos  0^  - (L  + r)j  1, 


(42) 


it  is  permissible  without  loss  of  accuracy,  ;o  let 


R ~ L + r1 
s 


(43) 


Then,  with  the  substitution 


sin  0 d0 
R 


dR 
r r1 


(44) 


the  integral  may  be  written 


Iw  = e 
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(45) 


where  0 is  regarded  as  a function  of  R. 


Integration  by  parts  now  gives 
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dR 


The  second  term  could  be  repeatedly  integrated  by  par+s  but  would  only  pro- 
duce terms  of  higher  order  in  k±  * . Consequently,  it  is  neglected  and  since 
the  first  term  is  zero  at  the  upper  limit,  the  integral  reduces  to 


±iks(L+r')  _a  (r.rl)  ik±(r-r') 


r r' 


ik, 


(47) 


The  above  approximation  is  valid  for 


k±  r<  k± 


k 

k±  o 


The  x integration  in  the  nearfield  is  represented  by 
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(49) 


With  the  approximation 


Rsl^  L + r'  ' 


the  integral  may  be  evaluated  analytically. 

The  variable  of  integration  is  transformed  to  R^  by  the  relation 
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With  a treated  as  small  compared  to  k±,  the  result  becomes 
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The  above  expressions  used  to  reduce  the  integration  to  a single 
range  integration,  have  produced  excellent  results.  Because  of  the  con- 
vergence problems  already  noted,  the  single  integration  program  is  believed 
to  be  inherently  more  accurate  than  the  double  integration  version,  at  least 
for  pump  frequencies  much  higher  than  the  signal  frequency. 


The  scattering  integral  and  one  dimensional  models  give  nearly 
the  same  results,  demonstrating  that  the  handling  of  the  finite  amplitude 
taper  is  consistent.  The  scattering  integral  approach  is,  of  course,  more 
accurate  in  handling  geometry.  With  the  reduction  to  a single  integration 
the  scattering  integral  computer  program  is  very  fast  and  inexpensive  to  use. 

1 

Results  obtained  using  the  single  integration  program  are  compared  with 
available  experimental  results  in  the  following  section. 

1. 

2.  3 Comparison  of  Theoretical  Results  with  Experimental 

Data  of  Other  Investigators 

Only  two  sources  of  experimental  data  showing  the  effects  of 

pump  saturation  on  parametric  receiver  axial  response  have  been  identified. 

8 

One  source  is  an  article  published  by  Berktay  and  Al-Temimi  which  contains 
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secondary  amplitude  response  curves  for  a pump  frequency  of  2850  kHz  and 
severed  different  signal  frequencies.  In  this  experiment,  the  receiving  hydro- 
phone was  placed  at  the  approximate  nearfield  collimation  distance  of  the 
pump  transducer.  The  other  source  is  experimental  data  obtained  by  Konrad, 
Mellen  and  Moffett.  ^ 

Earlier  efforts  by  the  present  author  to  compare  theory  with  the 

data  reported  by  Konrad,  Mellen  and  Moffett  led  to  confusion  because  even 

the  low  amplitude  data  did  not  agree  with  any  available  theory.  After  the 

discrepancies  were  brought  to  the  attention  of  the  authors,  Moffett  concluded 

that  a better  estimate  of  the  hydrophone  sensitivity  could  be  made.  He  did 

32 

so  and  was  kind  enough  to  provide  revised  data  , not  only  for  the  originally 
reported  experiment  which  used  a 620  kHz  pump  and  a 44  kHz  signal,  but 
also  for  experiments  performed  at  the  same  time  using  two  additional  signal 
frequencies,  30  and  13.  5 kHz.  It  is  the  revised  data  that  are  discussed  here. 

Theoretical  calculations  were  made  using  the  scattering  integral 
approach  detailed  in  Section  2.2.2.  Finite  amplitude  effects  were  accounted 
for  by  use  of  the  taper  function  defined  in  Section  2.2.  1.  The  data  so 
obtained  are  represented  by  the  solid  curves  in  Figures  2 through  5. 

Figure  2 shows  the  comparison  with  a typical  data  set  from  Berktay 
and  Al-Temimi.  Figures  3 through  5 show  comparisons  with  data  of  Konrad, 
Mellen  and  Moffett  for  two  different  ranges,  5.3  m and  9.1m,  and  three 
different  signal  frequencies.  The  data  for  the  44  kHz  signal  frequency  re- 
place the  ones  originally  reported. 

Figure  2 shows  a definite  disagreement  between  theory  and 
experiment  for  the  data  of  Berktay  and  Al-Temimi,  although  the  theory  does 
fit  the  low  amplitude  data.  Barring  experimental  error  these  data  indicate 
that  the  present  taper  function  results  in  too  much  loss  due  to  saturation  for 
the  interaction  in  the  nearfield  of  the  pump.  The  data  actually  agree  closely 
with  the  theory  given  by  Bartram  which  gives  similar  results  to  what  is 
obtained  by  the  present  methods  if  the  taper  function  is  taken  to  be  equal 
to  the  fundamental  component  amplitude  of  the  pump  signal,  rather  than  using 
the  taper  function  derived  in  Section  2.2. 


Secondary  Received  Level  Secondary  Received  Level 

- dB  re  1 p.Pa  (rms)  - dB  re  1 p.Pa  (rms) 


Pump  Source  Level  - dB  re  1 M-Pa  at  lm  (rms) 


Pump  Source  Level  - dB  re  1 H-Pa  at  lm  (rms) 
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Figure  4.  Comparison  of  Theory  and  Experimental  Data  of  Konrad, 

Mellen  and  Moffett  (Ref  32).  Pump  Frequency  = 620  kHz, 
Signal  Frequency  = 30  kHz,  Source  to  Hydrophone  Range 
= 21.  5m,  Rq  = 3.  25m,  Signal  Level  at  Hydrophone 
= 147.  2 dB  re  1 ^Pa 


(a)  Pump  to  Hydrophone  Range  = 5.  3m 

(b)  Pump  to  Hydrophone  Range  = 9.  lm 


Pump  Source  Le~rel  - dB  re  1 pPa  at  lm  (rms) 


Pump  Source  Level  - dB  re  1 pPa  at  lm  (rms) 

• EXPERIMENTAL  DATA 

CALCULATED  - SCATTERING 

INTEGRAL  MODEL 

Figure  5.  Comparison  of  Theory  and  Experimental  Data  of  Konrad, 

Mellen  and  Moffett  (Ref  32).  Pump  Frequency  = 620  kHz, 
Signal  Frequency  = 13.5  kHz,  Source  to  Hydrophone 
Range  = 21.  5m,  R0  = 3.  25m,  Signal  Level  at  Hydro- 
phone = 149.  8 dB  re  1 [J-Pa 

(a)  Pump  to  Hydrophone  Range  = 5.  3m 

(b)  Pump  to  Hydrophone  Range  = 9.  lm 


In  contrast  to  the  Berktay,  AI-Temimi  results  the  Konrad,  Mellen 
and  Moffett  data  support  the  present  theory  quite  well.  Figure  3,  which 
gives  the  44  kHz  signal  frequency  data,  shows  excellent  agreement  at  the 
9.  lm  hydrophone  range  and  good  agreement  so  far  as  curve  shape  is  con- 
cerned, but  a 7 dB  offset  at  the  5.3  m hydrophone  range.  Because  the  offset 
does  not  appear  in  any  of  the  five  other  data  sets  it  is  undoubtedly  the  result 
of  some  bias  error  in  the  experiment.  Figure  4 shows  quite  good  agreement 
between  theory  and  experiment  for  the  30  kHz  signal  frequency  experiment. 

It  is  interesting  to  note  that  for  either  range,  if  the  44  kHz  data 
are  overlayed  on  the  30  kHz  data  and  allowed  to  move  vertically,  the  two 
sets  may  be  made  to  line  up  almost  perfectly.  That  is,  the  experimental 
saturation  curves  have  the  same  shape.  This  shape  represents  slightly  less 
saturation  than  the  theory,  particularly  at  the  5.3  m range,  but  the  agreement 
is  considered  quite  good.  Here,  particularly  at  the  9.  lm  range,  the  present 
theory  gives  a much  better  fit  than  one  that  uses  the  fundamental  component 
amplitude  as  the  taper. 

Figure  5 gives  the  data  for  a signal  frequency  of  13.  5 kHz.  These 
data  show  a shape  considerably  different  from  that  of  the  other  sets.  Since 
only  the  signal  frequency  was  changed, and  there  is  no  known  reason  for 
this  to  appreciably  affect  the  shape  of  the  saturation  curve,  it  is  believed 
that  the  difference  may  have  been  caused  by  the  electronics.  Therefore, 
considerably  more  weight  is  placed  on  the  excellent  fits  obtained  in  the 
44  kHz  and  30  kHz  data. 

It  is  concluded  that  the  data  of  Konrad,  Mellen  and  Moffett 
generally  support  the  theory  presented  in  Section  2. 2 but  that  the  Berktay 
and  Al-Temimi  data  do  not.  The  derivation  of  the  theory  is  based  on  the 
assumption  of  farfield  propagation  and  the  possibility  exists  that  the  near- 
field may,  in  general,  exhibit  a different  type  of  behavior.  This  may  account 
for  the  different  behavior  evidenced  in  the  Berktay,  Al-Temimi  data.  This 
possibility  will  be  discussed  further  in  Section  2.5. 
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2.4 


Experimental  Measurement  of  P rametric 
Receiver  Amplitude  Response 


Insofar  as  directivity  r ieasurements  are  not  required  in  the  study 
of  amplitude  response,  it  is  possible  to  set  up  a well  aligned  and  controlled 
experiment  in  a tank.  Such  an  experimental  arrangement  was  effected  in  the 
AMETEK  test  pool  and  is  shown  schematically  in  Figure  6.  A 3.  8 x 11.  4 cm 
pump  transducer  was  used,  operating  at  a frequency  of  700  kHz.  Typical 
beam  patterns  for  this  projector  are  shown  in  Figures  7 and  8.  A 50  kHz 
signal  source  was  positioned  at  one  end  of  the  pool  7.  9 m behind  the  700  kHz 
pump.  The  narrow  beam  pump  projector  was  driven  by  a 200  watt  power 
amplifier  from  Scientific  System  Technology,  Inc.  The  pump  signal  and 
intermodulation  products  at  650  and  750  kHz  were  received  by  an  E-8  hydro- 
phone at  a distance  of  3.7m  from  the  pump.  The  50  kHz  signal  was  quite 
intense  by  comparison  with  the  usual  situation  in  parametric  receivers. 
Consequently,  it  was  quite  easy  to  separate  the  sidebands  from  the  pump 
by  narrow  band  filtering.  The  pump  was  operated  in  a pulsed  mode,  a com- 
promise necessitated  by  the  transducer  design.  However,  the  precaution  of 
using  a long  enough  pulse  to  ensure  the  gate  sidebands  being  well  below 
the  secondary  signals  was  taken.  A 4 ms  pulse  proved  to  be  adequate  and 
was  used. 


The  narrow  band  filtering  was  accomplished  by  passing  the  com- 
posite signal  into  a Hewlett-Packard  Model  310A  wave  analyzer  set  on  a 
bandwidth  of  200  Hz.  The  output  signal  was  passed  to  standard  calibration 
equipment  and  measured  on  a logarithmic  plotter.  The  pump  and  the  sum 
and  difference  frequencies  were  measured  by  peaking  the  wave  analyzer 
response  on  each  frequency  in  turn.  The  measurements  were  made  absolute 
by  passing  a 4 ms  pulse  of  known  amplitude  at  700  kHz  through  the  signal  path. 


Only  one  set  of  data  were  obtained  before  the  transducer  failed 
due  to  being  overdriven  and  efforts  to  build  a second  unit  were  unsuccessful. 
The  experimental  data  are  shown  in  Figures  9 and  10.  Figure  9 shows  the 
pump  fundamental  component  amplitude  response.  The  solid  cuwe  is  theory 
(Ref.  24).  These  data  are  quite  consistent  and  show  a definite  nonlinearity 


Figure  6.  Experiment  Configuration  - Parametric  Receiver 
Amplitude  Response 


SUM  FREQ  THEORY  - 
SCATTERING, INTEGRAL  MODEL 


at  high  source  level  in  agreement  with  theory.  Figure  10  shows  the  secondary 
frequency  amplitude  data.  Although  these  data  are  not  as  consistent  as  those 
for  the  pump,  agreement  with  the  theory  is  indicated.  The  solid  curve  is 
based  on  the  one  dimensional  theory  (Section  2.1).  Also  shown  on  the  plot 
are  points  calculated  by  the  scattering  integral  method. 

In  the  experimental  measurements  the  source  level  of  the  50  kHz 
signal  was  190.  5 dB  re  liiPa  at  1 m.  This  was  determined  by  pulsing  the 
source  and  measuring  the  amplitude  prior  to  the  onset  of  reflections.  In 
the  actual  parametric  receiver  experiment  the  signal  was  continuous  and 
considerable  standing  waves  were  present.  In  spite  of  this  the  data  agree 
with  theory  based  on  the  free-field  source  level.  The  receiving  array  is 
insensitive  to  reflections  not  collinear  with  the  pump. 

2 . 5 Conclusions:  Parametric  Receiving  Array  Finite  Amplitude  Taper 

Two  models  for  predicting  parametric  receiving  array  axial 
response  have  been  derived.  One  is  based  on  a one  dimensional  weak  shock 
theory  model  with  a simple  correction  for  nearfield  geometry.  The  other  is 
based  on  the  scattering  integral  approach. 

Both  models  have  weak  shock  theory  as  the  basis  for  predicting 
saturation  effects.  The  main  result  of  the  weak  shock  theory  analysis  is 
that  the  pump  signal  is  distorted  by  finite  amplitude  effects.  At  the  same 
time  each  distortion  component  of  the  pump  signal  is  frequency  modulated 
by  the  signal  frequency.  As  the  pump  fundamental  component,  which  is  the 
only  one  resulting  in  sum  and  difference  frequencies,  attenuates,  the  side- 
bands are  also  attenuated.  Thus,  in  the  scattering  integral  approach,  the 
taper  function  is  not  a positive  definite  quantity  but  may  become  negative, 
representing  a diminuation  of  sideband  energy  due  to  pump  saturation. 

This  view  of  parametric  receiving  array  taper  seems  physically 
reasonable,  at  least  for  one  dimensional  propagating  waves.  It  is  supported 
experimentally  by  the  results  of  Konrad,  Mellen  and  Moffett  and  by  the 
experimental  results  obtained  during  the  present  investigation,  although  the 
pump  amplitudes  obtained  in  the  latter  work  were  not  high  enough  to  make 
the  results  really  definitive. 
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The  present  taper  theory  applied  to  the  pump  nearfield  gives 
results  that  sharply  disagree  with  data  presented  by  Berktay  and Al-Temimi. 
Because  application  of  the  theory  to  the  nearfield  is  not  so  readily  justifiable 
as  its  application  to  the  farfield,  this  disagreement  is  taken  as  an  indication 
that  the  theory  should  only  be  applied  to  cases  in  which  the  receiving  hydro- 
phone is  well  in  the  pump's  farfield.  The  data  that  do  support  the  theory  are 
from  experiments  ir’  which  the  receiving  hydrophone  was  well  in  the  farfield 
of  the  pump. 


The  data  of  Berktay  and  Al-Temimi  tend  to  support  the  hypothesis 
that  the  taper  in  the  nearfield  should  be  the  amplitude  of  the  pump  fundamen- 
tal component  as  a function  of  range.  It  is  tempting  to  speculate  that  this 
may,  in  fact,  be  the  nearfield  taper  and  that  this  should  be  combined  with 
the  present  farfield  theory  to  give  a more  comprehensive  model.  The  present 
author  has  refrained  from  doing  this  because  the  justification  would  be  purely 
empirical  and  there  is  only  one  available  supporting  experiment.  Should  any 
additional  nearfield  data  become  available  the  question  of  what  taper  applies 
to  the  nearfield  should  again  be  addressed.  In  the  meantime  the  present 
model  is  believed  to  adequately  describe  the  axial  response  of  parametric 
receivers  in  which  the  receiving  hydrophone  is  well  in  the  farfield  of  the  pump. 
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3 . 0 TAPER  FUNCTIONS  AND  THE  PARAMETRIC  TRANSMITTING  ARRAY 

The  scattering  integral  approach  has  been  the  principal  method 
of  predicting  parametric  transmitting  array  performance,  whether  by  numerical 
volume  integration  or  by  approximating  the  integral  by  come  simpler  form.  In 
this  approach  the  concept  of  a taper  function  is  useful  because  it  enables  the 
finite  amplitude  effects  to  be  described  in  terms  of  a range  and  amplitude 
dependent  diminuation  of  the  source  strength  density. 

Taper  functions  have  previously  been  applied  successfully ^only 
to  simplified  integration  models  - for  example,  the  Mellen  and  Moffett  model. 

It  was  the  objective  of  the  present  work,  to  apply  the  taper  concept  to  the  full 
volume  integration  approach.  It  would  then  be  seen  whether, by  more  detailed 
account  of  geometry  and  inclusion  of  the  directivity  factor  in  calculating  the 
taper  function,  the  results  woulddiffer  significantly  among  various  taper 
functions  and  whether  they  would  differ  from  the  results  obtained  by  Mellen 

and  Moffett. 

Several  taper  functions  have  been  compared,  including  a new  one 
derived  using  weak  shock  theory.  It  should  be  noted,  however,  that  none  of 
the  taper  functions  considered  includes  account  of  possible  higher  order  inter- 
actions that  may  cause  the  secondary  frequency  to  be  depleted  at  extremely 
high  amplitudes.  Such  effects  are  included  in  a one  dimensional  model  pro- 
posed by  Fenlon.17  That  model  is  compared  with  the  results  of  Mellen  and 
Moffett  in  Ref.  13  wherein  it  is  shown  that  significant  differences  are  expected 
at  extremely  high  amplitudes.  In  order  to  include  these  higher  order  inter 
actions  in  the  scattering  integral  approach  it  would  be  desirable  to  derive  a 
taper  function  from  the  Fenlon  results  in  a manner  similar  to  the  treatment  of 
the  parametric  receiving  array  in  the  previous  sections.  However,  that  is 
beyond  the  scope  of  the  present  work. 

In  Section  3.  1 the  derivation  is  given  for  a taper  function  based 
on  a time  domain  theory  for  difference  frequency  generation  given  by  Merklinger . 
Because  the  taper  is  based  on  the  primary  intensity  it  will  be  referred  to  as  the 


I 

"Intensity  Taper  Function.  " As  mentioned  in  Section  1.  0,  the  phase  reversed 
parametric  array  data  of  Mellenbruch  and  Muir^re  useful  for  trying  out  a 
taper  function  concept,  Therefore,  the  phase  reversed  version  of  the  taper 
function  is  also  derived. 


* ' 


n 


1-4 


Numerical  results  are  compared  with  the  phase  reversed  data  in 
Section  3.2.  In  Section  3.  3 several  of  the  proposed  taper  functions  are  com- 
pared Results  of  parametric  array  calculations  are  compared  in  Section  3.  4. 


3 . l Derivation  of  the  Intensity  Taper  Function 

Merklinger^has  shown  that  the  source  function  for  secondary 
radiation  in  the  lower  sideband  is  proportional  to  the  time  derivative  of  the 
total  intensity  under  the  modulation  envelope.  The  difference  frequency 
source  function  is  obtained  by  extracting  the  Fourier  coefficient  at  the  dif- 
ference frequency 
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where 


(55) 


The  intensity  variation  in  the  envelope  is  assumed  to  be  unaffected 
by  the  presence  of  the  lower  sideband.  Therefore,  if  the  primary  spectrum 
consists  of  two  frequencies  of  equal  amplitude,  I(y)  may  be  calculated  from 
theory  for  an  initially  monochromatic  wave  at  the  carrier  frequency  of 
intensity  amplitude  IQ(y).  This  is  done  using  a modified  form  of  weak  shock 
theory  described  by  Lockwood2 4in  which  the  ordinary  weak  shock  theory 
solution  for  Bj  (a)  is  used,  with  a for  directive  spherical  waves  replaced  by 
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and  the  fundamental  component  amplitude  B^(ct)  multiplied  by  e 
The  problem  of  obtaining  I(y)  for  the  ordinary  and  the  phase  reversed  para- 
metric arrays  is  now  reduced  to  finding  appropriate  weak  shock  theory 
expressions  for  the  intensity  variation  as  a function  of  a . 

The  intensity  in  an  initially  sinusoidal  wave  without  reflection 
is  obtained  by  assuming  that  continuous  parts  of  the  wave  are  described  by 


V = sin  where 

0 = u (t-x/c)  + a sin  0 = y + ct  sin  0 . 

c 

The  dependence  of  the  total  intensity  is  then 
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It  is  convenient  to  normalize  the  preceding  expression  to  a zero  distortion 
value  of  unity.  Then  we  have 
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An  expression  for  the  intensity  in  an  initially  sinusoidal  wave 
with  reflection  at  a range  corresponding  to  j = a is  now  derived.  It  is 
assumed  that  after  inversion  straight  lines  remain  straight.  Otherwise, 
continuous  parts  are  described  by 

V = sin  0 , where 

0 = y + o-£  sin  0 , and  (61) 


We  redefine  0 , as  the  shock  phase  corresponding  to  the  shock  that  forms 
min 

subsequent  to  reflection;  0 is  the  phase  corresponding  to  the  pre-reflection 

max 

shock  0 =7 r - 0 . . (Note  that  the  time  origin  has  been  shifted  by  * 

' max  min  r 

to  put  the  new  shock  in  the  usual  place.) 

We  now  have  the  intensity  given  by 
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The  normalized  expression  in  terms  of  0max  and  0min  is 
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In  case  0 5 0 , the  preceding  expression  breaks  down.  The  straight 

*max  min 

segment  corresponding  to  the  pre-reflection  shock  overtakes  the  zero  and 
a sawtooth  is  formed  instantly.  The  normalized  intensity  turns  out  to  be 
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The  present  taper  function  (Eq.  60)  is  compared  with  other  proposed 
taper  functions  in  Section  3.  3.  First,  however,  results  computed  using  the 
taper  function  (Eqs.  63  and  64)  are  compared  with  experimental  data  for  the 
phase  reversed  parametric  array. 

3.  2 Numerical  Results  - Phase  Reversed  Parametric  Array 

7 

The  numerical  integration  program  of  Muir  and  Willett  has  been 
cast  into  a form  that  accepts  a finite  amplitude  taper  function.  The  taper 
function  which  is,  of  course,  amplitude  dependent,  is  calculated  as  a 
function  of  the  angular  coordinate  so  that  the  directivity  pattern  of  the 
primary  radiation  weights  the  amplitude.  Thus,  the  axial  taper  is  much 
more  extreme  than  the  taper  at,  for  example,  the  -3  dB  points. 


Several  methods  of  accounting  for  nearfield  saturation  effects 
have  been  tried.  The  best  results  have  been  obtained  by  treating  the  near- 
field as  a collimated  plane  wave  out  to  the  Rayleigh  distance.  The  ex- 
pression for  o is  given  in  Eq.  (70)  of  Section  3.  3.  It  is  interesting  to  note 
that  if  spherical  wave  theories  are  used  the  same  nearfield  distortion 
(assuming  aRQ  « 1)  may  be  obtained  by  letting  rQ  = RQ/e  . 
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As  a preliminary  calculation  to  check  out  the  program  and  the 

taper  function  derived  in  Section  3.  1,  theoretical  curves  were  calculated 

to  compare  with  the  phase  reversed  parametric  array  data  of  Mellenbruch 

and  Muir\  The  results  are  shown  in  Figure  11.  The  data  are  for  a 14  kHz 

difference  frequency  produced  by  a source  operating  at  424  and  438  kHz 

and  reflected  at  a range  of  4 yards.  Results  using  a less  sophisticated 

2 8 

theory  were  previously  reported  by  Lockwood.  The  present  results  agree 
considerably  better  with  experiment  than  the  previous  ones  did.  However, 
the  agreement  at  the  highest  amplitudes  is  not  as  good  as  had  been  hoped 
for.  Anomalous  nearfield  effects  are  suspected  as  the  reason. 

3.  3 Comparison  of  Various  Taper  Functions 

17 

Fenlon  has  given  a list  in  consistent  notation  of  several  taper 
functions  that  have  been  proposed  by  various  authors  for  use  in  conjunction 
with  quasi  one  dimensional  geometrical  models.  Those  from  Fenlon's  list 
that  apply  to  the  farfield  of  the  parametric  source  are  now  compared. 

First,  it  is  noted  that  all  of  the  functions  in  the  list  may  be 
expressed  in  the  form 
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(65) 


provided  that  u be  appropriately  defined.  The  form  of  Eq.  (65)  was  proposed 
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independently  by  Mellen  and  Moffett  and  by  Merklinger  as  an  approxi- 

18 

mation  of  the  fundamental  component  taper  given  by  Blackstock.  Equation 
(65)  is  compared  with  the  Blackstock  result  in  Figure  12,  which  also  shows  , 
for  a comparable  definition  of  cr,  the  taper  function  derived  in  Section  3.1. 
From  Figure  12  it  appears  that  the  differences  among  the  three  functions  are 
not  great. 


The  expressions  for  a appropriate  for  the  various  authors'  treat- 

13 

ments  of  Eq.  (65)  are  now  discussed.  The  Mellen- Moffett  form  is  given  by 


crQ  sinh"1  (r/rQ) 


(66) 
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Figure  11, 


Reflected  Parametric  Transmitter.  Reflector  at  4 yd. 
Difference  frequency  level  data  from  Mellenbruch  & 
Muir  (Ref.  1). 
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The  Merklinger,  Mellen  and  Moffett  taper  function,  hereafter  referred  to 
for  brevity  by  the  name  of  the  first  author  only,  is  given  by  Eq.  (65)  with 
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Two  taper  functions  due  to  Fenlon  are  the  Fenlon  Taper  Function  and  the 
Modified  Berktay,  Leahy  Taper.  The  appropriate  forms  of  a are,  for  Fenlon, 
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and  for  Modified  Berktay,  Leahy, 
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Here,  E^(-)  represents  the  exponential  integral.  Equation  (69)  is  equivalent 
to  the  ct  used  in  the  taper  function  derived  in  Section  3.  1 (Eq.  56). 


The  results  of  evaluating  Eq.  (65)  as  a function  of  r/rQ  with  a 
given  variously  by  Eqs.  (66),  (67),  (68),  and  (69)  are  shown  in  Figure  13. 
This  comparison  was  made  for  a case  when  atQ  rQ  was  small  (a0r0  = .001). 
Also,  for  the  Mellen-Moffett  calculation  the  value  of  rQ  was  taken  to  be 
twice  the  value  used  in  the  others.  It  will  be  seen  that  there  is  very  little 
difference  between  the  results  of  Eqs.  (67)  and  (68).  The  results  of  Eq.  (69) 
are  so  nearly  identical  to  those  of  Eq.  (68)  that  they  have  not  been  shown. 
The  results  using  Eq.  (66)  are  nearly  the  same  as  the  others  except  for  small 
and  large  values  of  r/rQ  . 

For  larger  values  of  Q r,  the  disparity  among  the  various  forms 
increases.  However,  it  is  likely  that  the  predominant  taper  will  be  absorp- 
tion in  cases  when  there  is  a significant  difference. 
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Figure  13.  Taper  Function  Comparison. 


3.  4 Comparison  of  Parametric  Transmitting  Array 

Calculations  Using  Various  Taper  Functions 

In  the  previous  section  the  relative  amplitudes  of  several  taper 

functions  were  compared.  The  taper  functions  considered  were  identified 

13  16  16  33 

as  the  Mellen-Moffett  , Fenlon  , modified  Berktay-Leahy  , Merklinger,' 

and  intensity  taper  functions.  It  was  noted  that  all  except  the  last  may  be 

represented  by  the  same  function  or  a single  parameter  called  a.  It  is  the 

definition  of  a that  distinguishes  the  various  tapers.  The  values  of  the 

taper  function  in  the  farfield  were  compared  for  a particular  set  of  parameters 

for  various  theories  excluding  the  intensity  taper,  and  it  was  concluded 

that  the  differences  among  them  were  quite  small. 

The  results  of  using  the  various  tapers  in  the  same  numerical 
volume  integration  program  are  now  compared.  As  a basis  for  comparison 
we  have  taken  the  parametric  scaling  curves  given  by  Mellen  and  Moffett. 
These  curves  are  the  solid  lines  in  Figure  14.  A single  value,  50,  of  the 
downshift  ratio  and  a wide  range  of  values  for  scaled  source  level  and  for 
(2a  Ro)-1  were  chosen.  Each  taper  function  except  the  one  derived  in  the 
present  investigation  was  defined  in  two  ways;  with  and  without  directivity 
taken  into  account  in  defining  off-axis  values  of  cr. 


The  Merklinger  and  Fenlon  ct's  are  as  defined  in  the  previous 
report  with  directivity  included  by  replacing  by  ctqD(8).  The  modified 
Berktay-Leahy  a is  the  same  as  the  expression  used  in  the  intensity  taper 
function,  which  is 
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The  version  without  directivity,  of  course,  lias  D(0)  - 1. 
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As  an  initial  comparison  the  low  amplitude  data  for  (2^  Rq)  - 1, 
10,  100  and  1000,  are  plotted  as  squares  (Q)  in  Figure  14.  For  the  source 
level  utilized  all  of  the  taper  functions  give  the  same  results.  It  is  interest- 
ing that  the  numerical  integration  gives  virtually  identical  results  to  the 
Mellen-  Moffett  model  at  low  amplitude,  with  the  exception  of  the  point 
corresponding  to  U^)-1  = 1,  which  falls  about  3 dB  above  the  solid  curve. 

The  higher  amplitude  results  for  the  modified  Berktay-Leahy  and 
the  intensity  taper  functions  are  also  shown  in  Figure  14.  This  comparison 
demonstrates  the  difference  in  results  between  the  two  functional  forms  of 
the  taper  theories,  with  the  same  cr.  A difference  of  about  3 dB  is  apparent 
in  each  of  the  cases  shown.  The  modified  Berktay-Leahy  results  without 
directivity  appear  about  1.  5 dB  lower  than  the  results  with  directivity, 
except  for  the  nearfield  limited  case.  The  results  discussed  above  are  all 
within  about  3 dB  of  the  Mellen-Moffett  results. 


A set  of  rather  surprising  results  was  obtained  by  using  the 
original  Mellen-Moffett  taper  function  in  the  numerical  integration.  The 
calculated  parametric  efficiencies  were  within  1 dB  of  the  original  model  ^ 
results.  These  data  are  not  shown  in  the  figures. 

The  same  data  points  were  run  using  the  Fenlon  and  the  Merklinger 
taper  functions.  Unfortunately,  the  results  using  these  tapers  are  not 
directly  comparable  with  those  discussed  above  because  they  are  b.  ed  on 
spherical  wave  theories  and  do  not  account  for  nearfielo  taper.  T 
results  could  have  been  made  to  fit  the  results  of  oti  ’ nodels  by  appropriate 
choice  of  the  effective  source  radius  rQ.  The  result1--  tud  indicate  that  the 
Fenlon  and  Merklinger  results  were  within  about  1 dB  of  each  other. 

It  may  be  concluded  that  the  parametric  transmitting  array 
scattering  integral  is  tolerant,  with  only  minor  differences  in  results,  of 
vast  differences  in  approach  and  in  degree  of  approximation  of  both  geometry 
and  of  the  taper  function.  With  all  of  the  results  lying  within  3 dB  of  each 
otter  it  would  be  quite  difficult  to  determine  by  direct  experiment  which  of 
the  approaches  gives  the  best  results.  As  was  noted  in  the  previous  section 





' ' 
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however,  an  adaptation  of  the  intensity  taper  function  to  the  reflected  para- 
metric array  produced  much  better  results  than  the  present  author  s earlier 
adaptation  of  the  Mellen-Moffett  theory. 

It  should  be  noted  that  all  of  the  data  used  in  the  above  compari- 
sons applied  to  range  points  well  in  the  farfield  of  the  interaction  volume. 

This  condition,  which  is  required  for  validity  of  the  Mellen-Moffett  model, 
is  not  necessary  in  the  volume  integration  approach, 

3.  5 Conclusions  - Parametric  Transmitting  Array  Finite  Amplitude  Effects 

The  results  reported  in  Section  3.  0 indicate  that  the  simple  geometry 
model  proposed  by  Mellen  and  Moffett  gives  results  almost  identical  to  those 
obtained  by  volume  integration,  when  the  same  taper  function  is  used.  They 
also  demonstrate  that  accounting  for  absorption  and  directivity  in  defining 
the  taper  function  associated  with  nonlinear  effects  makes  very  little  difference 
in  the  results.  On  the  other  hand,  the  manner  of  describing  nonlinear  effects 
in  the  nearfield  and,  in  spherical  wave  theories,  the  choice  of  the  effective 
source  radius  rQ,  make  substantial  differences  in  the  results. 

The  effects  of  higher  order  interact  ons  have  not  been  addressed 
in  the  present  work.  These  would  be  expected  to  be  important  only  at  extremely 
high  intensities. 
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RECOMMENDATIONS  FOR  FUTURE  work 


4.0 

The  results  of  the  present  investigation  point  to  a clear  require- 
ment  for  additional  experimental  work.  It  is  unfortunate  that  the  objectrves 
of  the  experimental  portion  of  this  work  were  not  fully  real^  because  of 

the  surprisingly  low  efficiency  of  the  pro  cs  that  were  built. 

objectives  should  be  pursued  further. 

A definite  requirement  exists  for  concentrating  on  nearfield 
saturation  effects,  both  in  the  parametric  receiving  array  and  In  the  para- 
metric transmitting  array.  In  the  case  of  the 

there  is  a fundamental  question  yet  to  be  answered  regarding  fo 
the  source  strength  taper  in  the  nearfield.  The  data  presently  «»*£“• 
suggest  that  the  nearfield  and  farfield  behavior  may  be  distinctly  ddfe  • 
The  theory  presented  herein  agrees  well  with  farfield  data  but  departs 
drastically  from  nearfield  data.  Additional  nearfield  data  are  required 
check  the  one  presently  available  source  and  to  form  the  basis  for  a Pos 
composite  model  that  would  include  both  nearfield  and  farfield 

In  addition  to  further  investigation  of  nearfield  effects,  addi- 
tional work  on  the  parametric  transmitting  array  should  address ^^iving 
a taper  function  that  accounts  for  higher  order  interactions.  Such  taper 
function  would  be  useful  because  of  the  flexibility  that  the  scattering  inte- 
gral approach  permits  in  accurately  describing  various  geometries. 
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